##Problem-01##
#R-Code#
z=read.csv("d:/s.csv")
z
attach(z)
data= data.frame(T1,T2,T3,T4,T5)
data
sum((T1-mean(T1))^2)

bigvector = c(T1,T2,T3,T4,T5)
bigvector
CF = (sum(bigvector))^2/length(bigvector)
CF

total.ss = sum(bigvector^2)-CF
total.ss

between.ss = (sum(T1)^2)/6 + (sum(T2)^2)/6 + (sum(T3)^2)/6 +(sum(T4)^2)/6 +(sum(T5)^2)/6 - CF
between.ss

within.ss = total.ss - between.ss
within.ss

df.total = length(bigvector)-1
df.total

df.between = ncol(data)-1
df.between

df.within = df.total-df.between
df.within

between.ms = between.ss/df.between
between.ms

within.ms = within.ss/df.within
within.ms

F.value = between.ms/within.ms
F.value

1-pf(F.value,2,21)

##C-Program##

#include<stdio.h>
#include<conio.h>
#include<math.h>
main()
{
    int a[6]={499,731,450,632,457,551};
    int b[6]={633,595,508,580,517,583};
    int c[6]={511,648,677,573,639,615};
    int d[6]={438,415,555,417,449,517};
    int e[6]={563,679,656,613,522,631};
    int ia=0,ib=0,ic=0,id=0,ie=0;
    int ias=0,ibs=0,ics=0,ids=0,ies=0;
    int g,yijs,ct,n=30,tss,trss,ri=6,ess,dft=4,dfe=25;
    float f,trms,erms,ftab=2.76;
    int ta=0,tb=0,tc=0,te=0,td=0;
    int tas=0,tbs=0,tcs=0,tes=0,tds=0;
    for(ia;ia<6;ia++)
{
    ta=ta+a[ia];

}
printf("\nT1=%d",ta);
for(ib;ib<6;ib++)
{
    tb=tb+b[ib];

}
printf("\nT2=%d",tb);
for(ic;ic<6;ic++)
{
    tc=tc+c[ic];

}
printf("\nT3=%d",tc);
for(id;id<6;id++)
{
    td=td+d[id];

}
printf("\nT4=%d",td);
for(ie;ie<6;ie++)
{
    te=te+e[ie];

}
printf("\nT5=%d",te);
g=ta+tb+tc+td+te;
printf("\nGrand Total=%d",g);
for(ias;ias<6;ias++)
{
    tas=tas+(a[ias]*a[ias]);
}
for(ibs;ibs<6;ibs++)
{
    tbs=tbs+(b[ibs]*b[ibs]);
}
for(ics;ics<6;ics++)
{
    tcs=tcs+(c[ics]*c[ics]);

}
for(ids;ids<6;ids++)
{
     tds=tds+(d[ids]*d[ids]);
}
for(ies;ies<6;ies++)
{
    tes=tes+(e[ies]*e[ies]);
}

yijs=tas+tbs+tcs+tds+tes;
printf("\n\nSum of yij^2=%d",yijs);

ct=(g*g)/n;
printf("\n\nCorrection Term=%d",ct);
tss=yijs-ct;
printf("\n\nTotal Sum of Squares=%d",tss);
trss=(((ta*ta+tb*tb+tc*tc+td*td+te*te)/ri)-ct);
printf("\n\nTreatment Sum of Squares=%d",trss);
ess=tss-trss;
printf("\n\nError Sum of Squares=%d",ess);
trms=trss/dft;
printf("\n\nTreatment Mean Square=%f",trms);
erms=ess/dfe;
printf("\n\nError Mean Square=%f",erms);
f=trms/erms;
printf("\n\nF Calculated Value=%f",f);
printf("\n\nF Tabulated Value=%f",ftab);
if(f>ftab)
	printf("\n\nSo we may conclude that Treatment means significantly differ");
	else
	 printf("\n\nSo we may conclude that All treatments are equally effective");
}



##Excel##

##See the excel shhet##



##Problem-02
#R-Code
x=c(62,59,54,46,57,53,50,64,55,55,48,52,53,50,57,53,63,50,57,53,50,54,57,55,41,64,55,52,57,60,53,62,50,47,55,50,55,50,56,47,53,56,61,68,55,55,59,58)
x
n=length(x)
xbar=mean(x)
xbar
nboot=1000
tmpdata=sample(x,n*nboot, replace=TRUE)
tmpdata
bootstrap_sample = matrix(tmpdata, nrow=n, ncol=nboot)
bootstrap_sample
bsmeans = colMeans(bootstrap_sample)
bsmeans
deltastar = bsmeans - xbar
deltastar
d=quantile(deltastar,c(.05,.95))
d
CI = xbar - c(d[1],d[2])
CI
plot(CI)

##C-Program##

#include<conio.h>
#include <stdio.h>
#include <stdlib.h>
#include<math.h>
#define N 1000
#define M 48
 main()
{
FILE *pt;
float x[M],y[N][M],ss[N]={0},mean[N],dd;
float var[N],ss1[N]={0},p,temp,ll,ul;
int i, pick,n[M],j;
pt=fopen("boot21.txt","r");

for(i=0;i<M;++i)
{
fscanf(pt,"%f\n",&x[i]);
}
fclose(pt);
srand (time(NULL));

for(i= 0; i <N; i++)
{
for (j= 0; j < M; j++)
{
n[j]= rand() % M;
y[i][j]=x[n[j]];
ss[i]=ss[i]+y[i][j];
mean[i]=ss[i]/M;
}
}
for(i=0;i<N;i++)
{
printf("\nMean of the sample no %i be %.2f\n",i+1,mean[i]);
}
    for(i=0;i<N;i++)
    {
        for(j=i;j<N;j++)
        {
            if(mean[i] > mean[j])
            {
                temp=mean[i];
                mean[i]=mean[j];
                mean[j]=temp;
         }
        }
    }
    printf("\nSorted mean:");
    for(i=0;i<N;i++)
    {
     printf("Sorted mean   %.2f\n",mean[i]);
    }
printf("Input the value of significance level \n");
scanf("\n%f",&p);
 dd=round(N*(p/2));
 pick=dd;
ll=mean[pick];
ul=mean[ N-pick];
printf("Confidence Interval  = %.2f\t\t%.2f",ll,ul);
return 0;
getch ();
}


##Excel_Law##

=INDEX(sample,ROWS(sample)*RAND()+1,COLUMNS(sample)*RAND()+1)



##Problem-03
#R-Code
placebo<-c(9243,9671,11792,13357,9055,6290,12412,18806)
oldpatch<-c(17649,12013,19979,21816,13850,9806,17208,29044)
newpatch<-c(16449,14614,17274,23798,12560,10157,16570,26325)
z <- oldpatch - placebo
z
y <- newpatch - oldpatch
y
ratio <- function(x, y) mean(x) / mean(y)
print(theta.hat <- ratio(y, z))
n <- length(y)
nboot <- 1000
theta.star <- double(nboot)
for (i in 1:nboot) 
{
    k.star <- sample(n, replace = TRUE)
    z.star <- z[k.star]
    y.star <- y[k.star]
    theta.star[i] <- ratio(y.star, z.star)
}
mean(theta.star)
hist(theta.star)
abline(v = theta.hat, lty = 2)
# bootstrap estimate of bias
print(bias.hat <- mean(theta.star) - theta.hat)
# bias corrected estimate
theta.hat - bias.hat

##Excel##
#See the 2nd problem#



##Problem-04
#R-Code

set.seed(type_your_roll_last_two_digit)
n=1000
m1=20
x=rpois(n,m1)
x
mle=mean(x)
mle
v=var(x)
v
var1=var(x)*(n-1)/n
var1
xb=0
var2=0
for (i in 1:n)
{
xb[i]=(sum(x)-x[i])/(n-1)
var2[i]=(sum((x-xb[i])^2)-(x[i]-xb[i])^2)/(n-1)
}
m=mean(xb)
m
v=mean(var2)
v
j=n*mle-(n-1)*m
j
jv=n*var1-(n-1)*v
jv
ulim=j+1.96*(sqrt(jv))
ulim
llim=j-1.96*(sqrt(jv))
llim


